Mon. Not. R. Astron. Soc. OOP. [TlfTil C2006') Printed 5 February 2008 (MN I*TeX style file v2.2) 



A Possible Origin of the Mass— Metallicity Relation of 
Galaxies 



J. Koppen^'^'^*, C. Weidner^'^f and P. Kroupa^'^^ 

^ UMR 7550, Observatoire Astronomique, 11 rue de I'Universite, F-67000 Strasbourg, France 
^ Institut fiir Astronomie und Astrophysik der Universitdt Kiel, D-24098 Kiel, Germany 
^ International Space University, Pare d'Innovation, F-67400 Illkirch, France 

Departemento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, Av. Vicuna MacKenna 4860, 
Macul, Santiago, Chile 

^ Argelander Institut fiir Astronomie der Universitdt Bonn, Auf dern Hiigel 71, D-53121 Bonn, Germany 
® The Rhine-Stellar-Dynamical Network 



Accepted 5 February 2008. Received 5 February 2008; in original form 5 February 2008 



ABSTRACT 

Observations show that galaxies follow a mass-metallicity relation over a wide range of 
masses. One currently favoured explanation is that less massive galaxies are less able 
to retain the gas and stellar ejecta and thus may lose the freshly produced metals in the 
form of galactic outflows. Galaxies with a low current star formation rate have been 
found to contain star clusters up to a lower mass limit. Since stars are predominately 
born in clusters, and less massive clusters have been found to be less likely to contain 
very massive stars, this implies that in environments or at times of low star formation, 
the stellar initial mass function does not extend to as high masses as during high star 
formation epochs. It is found that the oxygen yield is reduced by a factor of thirty 
when the star formation rate is decreased by 3 to 4 orders of magnitude. With this 
concept, chemical evolution models for galaxies of a range of masses are computed and 
shown to provide an e xcellent fit to the mass-metallicity relation derived recently by 
iTremonti et al.l ([2004') . Furthermore, the models match the relation between galaxy 
mass and effective yield. Thus, the scenario of a variable integrated stellar initial mass 
function, which is based on the concept of formation of stars in clusters, may offer 
an attractive alternative or partial explanation of the mass-metallicity relation in 
galaxies. 

Key words: stars: mass function - ISM: abundances - ISM: evolution - galaxies: 
abundances - galaxies: evolution - galaxies: ISM 



1 INTRODUCTION 

It has been known for a long time that there exist fairly 
well-defined relations between the me tallicity and mass (or 
luminosity) among gas-rich galaxies l|Lequeux et al.l 1 19791 : 
Lskillman ct al. 1989) as well as among gas-poor ga.laxies 
ijFabeiil 19731 : iBrodie fc Huchralll99ll ) . IZaritskv et all (|l994l ) 
point out the similarity of these relations, which is suggestive 
of a common physical origin of this phenomenon. Recently, 
the we alth of data from the Sloan Digital Sky Survey per- 
mitted ITremonti et al.l (|2004l ) to study the mass-metallicity 
relation over a wide range of masses and metallicities for 
about 53000 star-forming galaxies. 
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ITremonti et al.l (|2004l ) interpret their data with a sce- 
nario of gas outflows which are stronger in dwarf galaxies. 
More massive galaxies are able to retain the gas longer and 
more effectively than low mass objects, thus the metallic- 
ity can build up to higher values because following genera- 
tions of stars will be formed in an enriched gas environment, 
while the low-mass objects lose their metals through galac- 
tic winds. Different amounts of dark matter can also assist 
in this process (tPekel fc Silk . 19861 ). iDe Lucia et~al] (|2004l ) 
show that the chemical evolution in hierarchical galaxy 
formation models are able to explain the observed mass- 
metallicity relation and other relations like the TuUy-Fisher 
relation. This is achieved by including feedback processes 
like galactic winds into their cosmological simulations. The 
formulation of the feedback includes free parameters such as 
the feedback efficiency or the yield. 

While galactic outflows are observed and this scenario 
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provides a satisfactory match of the observed relation, one 
should not yet exclude other explanations or contributions 
from other processes. This is particularly important be- 
cause the current galactic wind models still need to be 
individually-fitted to the galaxies and are not easily ex- 
plained by phy sical parameters of t he galaxy. Further more, 
recent works bv lLee et all l|2006l ) and lDalcantoiil l|2006l ) show 
that simple outflows of ISM gas have di f ficult ies reproducing 
the yields in dwarf galaxies. iLee et al] (|2006l ) used observa- 
tions with Spitzer to derive the mass-metallicity relation 
from a number of dwarf irregular galaxies. They conclude 
that the dispersion in the relation and the large variations 
in the effective yields "are difficult to understand if galactic 
superwinds or outflows are res ponsible for low m etallicities 
at low mass or luminosity." In lDalcantonI l|2006l ) a series of 
closed-box chemical evolution models including infall and 
outflow is presented. There it is shown that neither simple 
infall nor outflow can reproduce the observed low effective 
yields in low mass galaxies, but metal-enriched outflows can 
do. That is, effectively the freshly synthesised elements need 
to be removed from the matter cycle, which is, in principle, 
equivalent to reducing the number of massive stars. 



In this work we like to explore further the implica- 
tions of_thfi_^£ttlM±staxs_jorm_g^ in clus- 
ters (|Lada fc Ladall2003h . lKroupa fc Weidned (|2003l ) showed 
that the integrated galactic IMF (IGIMF) of all stars in 
a galaxy depends on the embedded cluster mass function 
(ECMF), in the sense that more massive clusters are likely 
to have more massive stars, because of their greater amount 
of available gas to for m stars. Based on the relation found by 
IWeidner et al.l (|2004l ') between the maximum mass of clus- 
ters in a galaxy and the current star formation rate (SFR) , 
IWeidner fc Kroupal ()2005l . l2006l ) show that the high mass 
part of the resulting IGIMF is strongly dependent on the 
star formation history (SFH) in a galaxy, being in general 
steeper than the IMF. This giv es an important influenc e 
on the production r ate o f metals l|Weidner fc Kroupall200(i ). 
iGoodwin fc Paeell (|2005h explored the effect of the IGIMF 
on the metal production in galaxies. They found a change 
in the general metal yield by a factor of ~ 1.8 and by a fac- 
tor of ~ 1.5 for oxygen and magnes ium when they use the 
ECMF slope of iLada fc Ladal (|2003l . i.e. /3 = 2). 



One thus expects that the effective metal yields from the 
average stellar population may vary with time and galaxy 
type. It is not the aim of this contribution to exclude the sce- 
nario of galactic outflows, but to explore the possibilities and 
consequences arising purely from clustered star formation, 
and to probe to what extend it may provide an alternative 
or supporting explanation. 

Section [2] briefly reviews the model of the SFR- 
dependent IGIMF (henceforth IGIMF for simplicity), sec- 
tion [3] discusses the importance of the oxygen yields for this 
study, while in section |4] the chemical evolution model is in- 
troduced and the effect of a variable IGIMF on the chemical 
evolution of galaxies is studied. In section[5]the implications 
of these results are discussed, while section [6] presents the 
conclusions from these results. 



2 THE VARIABLE IMF 



InlKroupa fc Weidneij (|2003h and I Weidner fc Kroupal (|2005l . 
|2006| ) it has been shown that the formation of stars pre- 
dominantly in clusters has a profound influence on the in- 
tegrated galactic stellar initial mass function (IGIMF) of 
all stars in a galaxy. In contrast to the IGIMF, which is 
the "IMF" integrated over a whole galaxy, the IMF of 
stars formed in a single cluster is well described by the 
well-known power -law function with the Salpeter/Massey 
slope a = 2.35 l|Salpeteij 1 19551 : iMassevI 12003) for stars 
more massive than 1 Mp,, as found in several observational 



studies (iMassev fc Hunted Il99 8': 'Siria nni et"al]l200d. 
Parker et al.l l200ll: lMassevl"l20 02. 20031: IWvse et al. " 
BeU et al.ll2003l : iPiskunov et al.ll2004l) . with 



2002 



2002 



^(m) oc m 



(1) 



where £^{m) dm is the number of stars in the mass inter- 
val m, m + dm. Below . 5 Mq the IMF is found to flat- 
ten (iMiller fc Scald Il979l : IScalol 1 19861 : iKroupa et al] Il993l : 
iReid et al.ll2002l ). and t he standard IMF, becomes a seg- 
mented power-law (|Kroup a 2001, 2002): 



ao = 0.30 
ai = 1.30 
aa = 2.35 



0.01 s; m/Me < 0.08, 
0.08 «; m/Me < 0.50, 
0.50 < m/M0, < 
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Note that this form constitutes a two-part power-law stel- 
lar IMF with the "Salpeter/Massey" (or short "Salpeter") 
slo pe above .5 M p,. Note also that this form is similar 
to iKennicuttI (|l983l ') two-part power-law form, except that 
the slope at the high-mass end is less-steep here (being the 
Salpeter/Massey value), and that the break in the power- 
law occurs at a smaller mass here. These differences, while 
fairly subtle, are important as the form used here is based 
on a very thorough star-couiit ana lysis using post-1983 data 
l|Kroupa et al.lll993l : lKroupall200ll '). The physical stellar up- 
per mass limit is assume d to be mmax* = 15 Mq, f ollow- 
ing th e recent results by Weidner fc Kroupal (|2004 ). iFigerl 
l|2005l ). IOev fc Clarke! (|2005l ) and lKoenI l|2006l ). 

The masses of young, embedded star clusters also fol- 
low a power-law distribution, the embedded cluster mass 
function (ECMF): 



(3) 



where ^oci(Meci) dM^ci is the number of embedded clusters 
in the mass i nterval (Med, M^d + dMgci), and Med is the 
mass in stars. iLada fc Ladal (|2003h flnd a slope /3 r; 2 in the 
solar neighbourhoo d for clusters with masses between 50 
and 1000 Mq, while iHunter et all (|2003) find 2</3<2.4 for 
lO^Meci/M0 < 10* in the SMC and LMC, IZhang fc Fail 
l| 19991 ) flnd 1.95±a03forl0* ^M^i/Mq < 10'^ in the Anten- 
nae galaxies. IWeidner et all (|2004| ) find that (3 = 2.35 best 
reproduces the observed correlation between the brightest 
young cluster and the global star formation rate for a sam- 
ple of late- type galaxies. 

The IMF of all stars ever formed in a galaxy (the 
IGIMF) is obtained by summing up all standard IMFs of 
all clusters from which a galaxy is formed, 

/•A4cl,max 

6GiMF(m) = / ^(m ^ (Med)) ^ocl(Mecl) dMed, (4) 
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Figure 1. The dotted line is the standard IMF $(m), given by 
Eq. [2] which has a = 2.35 for m > 0.5 Mq. The dashed line is 
the overall IGIMF, $igimf("^), for /3 = 2.2. All functions are 
scaled to have the same number of objects in the mass interval 
0.01 — 1.0 M0. Slopes a = 2.35 and 2.77 are indicated by lines. 
Note that we use the notation log = log^Q throughout the paper. 



log Mccl,max(t) = log /CML + 0.75 ■ log SFR(t) + 6.77, (5) 

and a mass-to-light ratio of /cml = 0.0144, typically for 
young (< 6 Myr) clusters, the IGIMF evaluated during the 
time interval t, t + 5t - with the (c onstant) "star formation 
epoch" St = 10 Myr (|Weidner fc K roupa 200i) - becomes: 

/■Mocl,max{SFR(t)) 

6GiMF(m,SFR(t)) = / ({m^ (Mecl))Cocl(Afecl)dAfecl 

A similar timescale St is found bv lEgusa et al.l l|2004l ) from 
measuring the offset of H ll regions from the dark spiral 
arms in a number of galaxies. In our interpretation, St is 
the time scale within which an ECMF is completely pop- 
ulated. Th i s resu lt is at least qualitatively supported by 
lElmegreenI ||2000D . By comparing the observational prop- 
erties of embedded clusters and molecular clouds he de- 
rives "that s tar formation oc c urs in only one or two cross- 
ing times". iHartmann et al.l l|200ll) also argue for short 
< few Myr molecul ar-cloud life times. And the Antennae 
IZhang fc Falllll999l ') have a fully developed cluster MF at 
age of ~ 10 Myr. Thus, a lower SFR results in a lower maxi- 
mum cluster mass which essentially reduces the upper mass 
limit for the stars. Note that eq.|6]is the correct "IMF" to be 
used when study ing the global SFR-properties of galaxies. 
iKennicutd (llQsd ). for example, finds a = 2.5 (for m > 1 Mq) 
reproduces integrated luminosities of star-forming galaxies. 
Although he refers to this "IMF" as a "Salpeter IMF" , it is 
actually a steeper function (with a larger a). 



where ^cci(Mcci) is the ECMF and ^(m ^ mmax(Afoci)) is the 
stellar IMF in a particular cluster within which the maximal 
mass of a star is mmax. A critical ingredient of this concept is 
the existence of the correlation between the maximum stellar 
mass mmax and the cluster mass Mod, which depends on the 
sequence with which star s of different masses are formed. 
IWeidner fc Kroupal l|2006l ) find that several young clusters 
are best described if one assumes that, starting at the lowest 
mass, stars of progressively higher mass are formed. 

In Fig. [1] we show an example of the resulting IGIMF 
from Eq. 21 using an ECMF power-law slope P = 2.2. 
Note that the IGIMF {dashed line), which has a slope 
above 1.5 of about 2.77, is steeper than the stan- 
dard IMF (dotted line); furthermore it turns down rather 
sharply near mma x* = 150 M q , the fundamental upper mass 
limit for stars (Weidner fc K roupa 2004, 2005; Figcr 200^, 
lOev fcClaxQl2005yKoenll2006l'l. 



ElmegreenI 1 20061 ) puts forward the view that the 
IGIMF equals the IMF, but we note that his discussion 
points to a few examples which actually do impose a variable 
IGIMF. Our formulation tests, like that of Elmegreen, the 
case of l3 — 2, but ours allows a variation of the maximum 
mass of the stars with the SFR, in contrast to Elmegreen's 
where this mass comes out to be constant. 



2.1 The connection to the star formation rate 

The IGIMF depends not only on th e ECMF but also on 
the star formation history of a galaxy. IWeidner et al.1 (|2004l ) 
derive a relation between the current SFR and the maximum 
mass of a cluster in a galaxy. With this relation 



2.2 Model Parameters 

Because the IGIMF is steeper than Salpeter's IMF, which is 
often used in chemical evolution studies, its dependence on 
the SFR can be expected to have observable consequences 
for the chemical enrichment of galaxies. For our evolution 
models we shall use the following parameters: 

• minimum mass of the clusters: Afcci.min = 5M0. This 
corresponds to a Taurus- Auriga star-forming "cluster" con- 
taining a dozen stars, 

• range of "stellar" masses: 0.01 to 150 Mq, 

• ECMF power-law exponent /3 = 2 and also 2.35, 

• the standard IMF (Eq.[2j. 

The IGIMFs obtained for various SFRs are displayed 
in Fig. [5] in the form of number of stars per unit stellar 
mass interval (Salpeter's IMF has the exponent 2.35). A 
lower SFR reduces the maximum cluster mass and hence 
the upper mass limit for the stars; this is quite a strong 
effect. One also notes the steepening of the slope, especially 
at high stellar masses. 

In Fig [3] we show that the slope /3 of the cluster mass 
function also affects the slope of the IGIMF above stellar 
masses of 1 Mq. For the same SFR, the upper mass limit 
does not change, however increasing (3 further steepens the 
IGIMF slope. To display this dependence, we took a larger 
range of values for l3 than indicated from observations. 



3 OXYGEN YIELDS 

From the IGIMFs we compute the yields of oxy- 
gen by integrating over the entire stellar mass spec- 
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Figure 2. The IGIMFs computed for SFR = 0.01, 0.1, 1, 10, 
and 100 Mq/jv (solid lines, from bottom to top) AND /3 = 2, 
compared to the Salpeter IMF for the entire stellar mass range 
0.01 to 150 Mq. All IMFs are normalised to have the same number 
of stars in the interval 0.01 to 150 Mq . The standard IMF is Eq.[2] 
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Figure 3. Like Fig. [2] but for the same SFR = 0.01 Mo/yr and 
P = 1.65,2.0,2.35, and 2.7 (solid lines, from right to left). 



trum t he stellar yields weight ed with the IGIMFs 
(cf. e.g. iKoppen fc Arimotol Il99ll ). The prescriptions for 
the stellar nucleo s ynthe sis are: for massive stars we use 
iThielemann et al] (1 19961). including th e yields for the 40 
Mq Star quoted bv lThomas et al.1 (| 19981 ): beyond 40 Mq we 
assume that the mass fractions of freshly produced ele- 
ments a re constant. For low and intermedia te mass stars 
we take Ivan den Hoek fc GroenewegenI (|l997l ). This choice 
is known to give a good match to the chemical evolution of 
the solar neighbourh ood, if one uses a Salpeter IMF (cf. e.g. 

omas et al.lll998l) , because the yield is close to the solar 
oxygen abundance. As will be shown by the results below 
(Fig. Ill|> this choice also gives a good match with the obser- 
vational data for galaxy masses of IO^^Mq. Uncertainties in 
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Figure 4. Oxygen yields as a function of the SFR for several 
values of /3. The yields for a Salpeter IMF over all stellar masses 
and the standard IMF (Eq. [2J are shown by horizontal dashed 
lines. 

the stellar yields w ould amount to at least a factor of 2 (e.g. 
iHenrv et al.ll2000l ) to be applied in either direction. 

The oxygen yield (Fig. U) decreases strongly when the 
SFR is below about 1 Mq /yr. This mirrors the strong depen- 
dence of the upper mass limit on the SFR. Test calculations 
show that the influence of the variation of the IGIMF slope 
with SFR is of less importance. Since for high values of the 
SFR the IGIMF approaches the standard IMF (cf. Fig. [2), 
the yield converges towards the value obtained with the full 
standard IMF sampled to the maximum cluster mass. The 
yields also depend on the assumed value for 13; larger val- 
ues result in IGIMFs more strongly cut-off at the high mass 
end. In the figure we show - to illustrate this dependence 
- a greater range of values for [3 than suggested by obser- 
vations (2.0 ... 2.35; cf. section [2]). A Salpeter IMF for the 
entire stellar mass range gives an oxygen yield very close to 
the solar oxygen abundance. As the oxygen yields from mas- 
sive stars depend on the stella r metallicity only very weakly 
(see lWooslev fc Weaveill 19951 ). our results computed for so- 
lar composition stars also apply to non-solar metallicities. 

Stellar oxygen yields of stars as massive as 40 MQare 
less certain due to the poor knowledge ab out the explo- 
sion characteristics jWooslev fc Weaver! 19951 ). Furthermore, 
our adopted nucleosynthesis recipe extrapolates the stellar 
yields beyond 70 Mq . We investigated the influence of these 
uncertainties by computing the yields under the assumption 
that oxygen is ejected only by stars below a certain upper 
mass limit. As depicted in Fig. [5] the resulting yields are re- 
duced by about a factor of two, if oxygen is produced only by 
stars less than 40 Mq . The dependence on the SFR becomes 
somewhat less steep for low values of the SFR. 

If all galaxies had the same age of about 10 Gyr, galax- 
ies with masses between lO'^ to 10^^ Mq have time-averaged 
SFRs ranging from 0.001 to 10 Mq/jx (cf. Fig.gl). Since over 
this interval, the oxygen yield increases by about a factor of 
thirty, we may thus expect that the SFR-dependent IGIMF 
causes observable effects in the chemical evolution of these 
systems. 
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Figure 5. Same as Fig. [J] but for = 2 and several values for 
the upper mass limit for oxygen production in massive stars (from 
150 to 25 Mq). 



4 CHEMICAL EVOLUTION MODELS FOR 
GALAXIES OF DIFFERENT MASSES 

In order to study the effects of the SFR-dependent IGIMF 
on the present-day gas abundances of galaxies of different 
masses, we compute a series of illustrative chemical evolu- 
tion models to display the principal effects and show their 
magnitude. These models use simple parameterisations, and 
thus are not designed to match the exact time histories of 
the galaxies, however we concentrate on those properties - 
metallicity and gas fraction - in which chemical evolution 
models are known to be less sensitive to modelling details. 
Thus, the results retain a certain degree of general validity. 



4.1 Physical Assumptions 

Each galaxy is treated as a one-zone model, dealing with 
the gas return and metal production of stars of different 
masses in detail and taking into account the finite lifetimes 
of the stars. The ejected gas and metals are assumed to 
be mixed perfectly and instantaneously with the gas in the 
model, as is done in the usual chemical evoluti on stud- 
ies (cf. e.g. iTalbot fc Arnett|[l97ll : lMatteuccilll986l ). As de- 
scrib ed before, the stellar nuc leosynthesis prescriptions are 
from lThielema,nn et al.l lll996l') for massive stars (typ e II su- 
pernovae) and van den Hoek &: GroenewegenI (Il997l ) for low 



and intermediate-mass stars. 

At each time-step the IGIMF is selected in accordance 
with the current value of the SFR, using a pre-computed 
lookup table. The SFR has the same form as used by e.g. 
iMatteucci fc Francois! l| 19921 ). 



SPJ^(^) ^ MgasW+ A/stars fc) 
TSFR 



A/gas (t) 



(7) 



which depends on the current masses of gas, Afgas(f), and 
stars, A^stars(i), and thus on the current values of total mass 
and gas fraction. The parameter x is taken as unity, although 
test calculations are done with other values, giving very simi- 
lar results, provided the present-day gas fraction remains the 



same, as one should expect for closed-box and infall mod- 
els. The star formation time scale tsfr is assumed to be the 
same for all galaxy masses; several values shall be taken, 
which will thus determine the present gas fraction. In later 
models, we shall allow tsfr to vary with galaxy mass. 

The galaxy is either treated as a closed box with all 
its initial mass in the form of metal-free gas, or as an infall 
model where the primordial gas is accreted with a rate that 
decreases exponentially in time: 



M{t) = Afgal 



exp(^t/rf) 
1 - exp(-T/r/) ' 



(8) 



with the time scale r f = 5 Gyr, similar to [ Chi osi (198CD, 
IMatteucci fc FrancoisI (|l992l ) and iThomas et all ((1998 ) for 
models of the solar neighbourhood. We assume all the mass 
of the gas to be accreted within the age T = 13 Gyr of 
the galaxy. It is worth pointing out that with this kind of 
monotonic gas infall, the oxygen abundance remains close 
to the values from a closed-box m odel of the same gas frac- 
tions l|K6ppen fc Edmunds! [19991 ) . Thus, the choice of the 
infall time scale is of minor importance in the context of the 
primary element oxygen ejected with very little time delay. 

For the numerical calculations we use a grid of 190 loga- 
rithmically spaced stellar masses and 500 time steps to cover 
13 Gyr. 



4.2 Closed-Box Models 

We first compute a series of closed-box models. The star 
formation time scale is taken to be a fraction or multiple of 
the age of the galaxy: tsfr = T/a. The factor a is assumed 
to take values of 0.5, 1, 2, 3, and 4, i.e. time scales of 26, 13, 
7, 3.5 and 1.75 Gyr. These five series of models with different 
galactic masses span a rather wide range of present-day gas 
fractions, 



/gas(T) = 



Mga.(r) 



Afgas(r) +M,tars(r)' 



(9) 



In Fig. |6] we com pare these models with the data of 
iTremonti etaU (|2004l ) for the relation between stellar mass 
of the galaxies and gas-phase oxygen abundance. The com- 
puted relations are somewhat flatter than the observed re- 
lation, but models with present-day gas fractions between 
25 and 50% are indeed well suitable to explain the observed 
metallicities, which correspond to a ~ 1. For comparison we 
show one model series done with a constant Salpeter IMF. 
Since in these models the oxygen abundance is a function 
of only the gas fraction which is determined by the star 
formation time scale, each of the five model-series give an 
oxygen abundance independent of the galactic mass. The 
model shown is done with o = 1, the other models result 
in horizontal lines at different oxygen abundances, ranging 
from 8.5 for a = 0.2 5 to 9 4 for a = 4. 

ITremonti et al.l (|2004h also determine the relation be- 
tween total galaxy mass Afgai = Afgas (T) + Afstars (T) and 
the effective oxygen yield, which they compute from current 
oxygen abundances and gas fractions by assuming a closed 
box model. 



In /gas (T) ■ 



(10) 
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Figure 6. The relation between oxygen abundance and stellar 
mass of galaxies. Dashed lines indicate the 16, 50, and 84 per- 
centile s of the observed metallicity distribution bv lTremonti et ahl 
The symbols are results from closed-box models: From 
top to bottom, the model series (with (5 = 2) correspond to the 
factor a in the SFR to 4, 3, 2, 1, and 0.5. The symbols indi- 
cate the present-day gas fraction: large circles > 0.5, small cir- 
cles 0.5. ..0.25, large filled circles 0.25. ..0.15, small filled circles 
0.15. ..0.1, and crosses < 0.1. The plus-signs show the model se- 
ries with a = 1 but /3 = 2.35. The full horizontal line refer to a 
model grid obtained with a constant Salpeter IMF between 0.01 
and 150 Mq. 



Our closed box model sequences (Fig. [T]) reproduce the slope 
of the relation remarkably well, although being higher than 
the median observed relation by 0.1 dex. This is in strong 
contrast to models with a constant IMF which all have the 
same effective yield for all galactic masses and all star for- 
mation time scales. 

4.3 Infall Models 

Motivated from chemical evolution models for the solar 
neighbou rhood, we consider continuo u s infall of prirnordial 
Sas (e.g. 



iMatteucci fc FrancoisI Il992l : Ixhomas etal] 1 19981 : 

lisfaab &: OstrikeT 20061) . The infall rate is assumed to de- 



crease exponentially in time with a time constant of 5 Gyr. 
As can be seen from Fig. [S] one obtains larger gas fractions 
and lower metallicities. Thus, models with a « 2 are quite 
suitable. We note that infall does not change the slope of 
the relation which remains slightly too flat. 

As expected for models with monotonically decreasing 
infall, the effective yields are only slightly (less than 0.1 dex) 
lower than found from the closed-box models. Hence, they 
also match remarkably well the shape of the observed me- 
dian relation. 

Since a general trend is observed that massive 
galaxies have smaller gas fractions than dwarfs (e.g. 
iKoppen fc Hensleijl2005l l. we construct another set of mod- 
els in which the star formation time scale is assumed to 
decrease with increasing galactic mass Mgai in such a way. 



TSFR = 20 Gyr 



' I lOSMg 



gal 



Figure 7. The relation between the effective O yields and the 
total galactic mass: dashed lines indicate the 16, 50, and 84 per - 
centiles of the observed distribution by iTremonti et al.l ll2004h . 
The symbols denote the same sequences of closed box models 
of Fig. \E\ The horizontal line refers to models with a constant 
Salpeter IMF. The plus-signs show the model series with a = 1 
but /3 = 2.35. 
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Figure 8. Same as Fig. |6] but for models with an exponentially 
decreasing infall of primordial gas. 



that the observed data (Fig. [9]) are reasonably well matched. 
In Fig. [TO] we show three model sequences, done with the 
above time scale, as well as with half and twice its value. 
For comparison, we also compute three model sequences 
with a constant Salpeter IMF. A similar dependence of 
the SFR time scale on galactic mass (tsfr = 50 Gyr • 



(Afgal/lO^Mo) 



gives a satisfactory match to the ob- 



(11) 



servations. 

The resulting mo del sequences ma t ch ve ry well the dis- 
tribution observed bv iTremonti et ahl (120041 ) . as is seen in 
Fig. 111! It is worth noting that the agreement is excellent 
both in terms of slope and metallicity. 

The effective yields (Fig. I12|) also match the observed 
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IOg(MH|/MsuN) log(MsTARs/MsuN) 



Figure 9. The obs erved relation of hydrogen mass and gas frac- 
tion (adapted from lK"oDpen &: Henslerll2005l . where details may 
be found). The filled dots mark objects with gas fractions de- 
termined by iGarnettI ll2002h that account for both atomic and 
molecular content. The dotted line is an assumed relation, used 
for Fig. [19] 



Figure 11. Similar to Fig. [6] but for models with exponentially 
decreasing infall of primordial gas and with a star formation rate 
depending on the galactic mass (Eqns.[7]and[TTJ. The three series 
of models are as in Fig. 1101 The full curves refer to a model com- 
puted with a constant Salpeter IMF, while the plus-signs show 
the model series with a = 1 but (3 = 2.35. 




log(MH|/MsuN) 

Figure 10. The relation between present gas mass and gas frac- 
tion from models with gas infall of 5 Gyr time scale. The central 
series of symbols refers to star-formation time scales which depend 
on the galaxy's mass following Eqn. 1111 The SFR time scales are 
halved in the upper series and doubled for the lower sequence. 
The meaning of the symbols is as in Fig. \E\ The three full lines 
are for similar models, but computed with a constant Salpeter 
IMF. 



relation very reasonably, although the slope among the 
massive galaxies is flatter and even turns down, and near 
IO^'-^Mq the yields are about 0.1 dex too large. 

The adopted infall time scale is not a critical parame- 
ter for this match with observations: monotonically decreas- 
ing infall leads to but minor deviations from the closed-box 
model properties, in particular the relation between gas frac- 
tion and metallicity. As is evident from Figs. [6] and |8] the 




log({MsTARs +Mgas )/Msun ) 

Figure 12. Similar to Fig. [T] but for the infall models of Fig. 
1101 The symbols are the same as in Fig. [6] The nearly horizontal 
curves are from the models computed with a constant Salpeter 
IMF. 

metallicity in more massive galaxies are decreased by about 
0.2 dex; but it is worth emphasising that the locus of the 
galaxies with the same gas fraction has hardly changed at 
all. The same holds true if one chooses the infall time scale 
to be anywhere between about 1 and 10 Gyr. 

The models with a constant IMF do not provide such 
a good match: While the assumed relation between SFR 
time scale and galactic mass could account for an increase 
of metallicity with galactic mass, our simple recipe cannot 
reproduce the shape of the observed relation (Fig. Ilip . Fur- 
thermore, these models fail completely to match the mass- 
effective-yield relation, as their yields are nearly indepen- 
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Figure 13. Like Fig. 1111 but assuming that for galaxies with a 
total mass below 10^" Mq star formation occurs in bursts of 100 
Myr duration every 1 Gyr. During a burst the SFR is 30 times 
higher than the average rate. 



Figure 14. Like Fig. 1121 but assuming that for galaxies with a 
total mass below 10^" Mq star formation occurs in bursts of 100 
Myr duration every 1 Gyr. During a burst the SFR, is 30 times 
higher than the average rate. 



dent of galactic mass (Fig. I12|l . Thus the dependence of ef- 
fective yield on galaxy mass might be more crucial in tracing 
the origins of the mass-met allicity relation. 

4.4 Models with starbursts 

Among dwarf irregular galaxies the star formation history 
is kn own to show strong temporal variations (e.g. iGrebell 
I2OOII ). How does the presence of star bursts influence the 
results? Let us consider the extreme case, that all star for- 
mation occurs in short periods with long intervals of no star 
formation. We shall assume that in galaxies with less than 
10^° Mo stars are formed in a burst lasting 100 Myrs every 
1 Gyr. In order to obtain the same present gas fraction as 
in a model with continuous star formation, the SFR during 
the short bursts must be 30 times higher than the average 
rate. The resulting IGIMF contains a larger fraction of mas- 
sive stars, and thus the average metal yield will be higher 
than under continuous star formation. Thus, strong star- 
bursts lead to both higher oxygen abundances (Fig. I13p and 
higher effective yields (Fig. I14|l. as shown in this somewhat 
extreme model. One notes that since for these models the 
SFR is higher, the variations of the upper stellar mass limit 
and the IGIMF slope should be milder. As Fig. [14] shows, 
the variation of the effective oxygen yield is less than in 
the models with continuous star formation (cf. Fig. I12p . In 
the oxygen abundances, this change is much less apparent, 
because the assumed dependence of the SFR time scale on 
galactic mass (Eqn. Ilip contributes fairly strongly to the 
mass-met allicity relation (cf. Fig. Ilip . 

We emphasise that this schematic modelling of star for- 
mation bursts serves to show the influence on the results of 
a rather extreme burst-type history of star formation, irre- 
spective of the physical processes which could form the base 
for such a behaviour. The parameters were chosen to pro- 
vide a clear demonstration of the effects. Nonetheless, the 
results predict that with a SFR dependent IGIMF one would 
expect that galaxies which had a strongly fluctuating star 



formation history would tend to have higher metallicities 
and higher oxygen yields. 

5 DISCUSSION 

5.1 Chemical evolution models 

The consequences of the IGIMF depending on the SFR is 
primarily that the oxygen yield is higher in environments 

- or times - where the SFR is high. Thus, more massive 
galaxies are expected to have higher yields and thus have a 
tendency to have higher metallicities. This basic implication 
is quantified in the chemical evolution models presented in 
Sect. H 

• The closed-box models exhibit mass-metallicity and 
mass-yield relations which are very similar to what is ob- 
served, without the need to invoke additional physics. 

• Adding infall of primordial gas modifies both relations 
only slightly. The presence of gas infall or the true history 
of the infall - as long as it decreases monotonically in time 

- are not crucial for this comparison. 

• If the observed trend between galaxy mass and present 
gas fraction is matched by adopting a dependence of the star 
fo rmation time sca l e on g alaxy mass, both observed relations 
of iTremonti et all l|20o4 i can be reproduced very well. This 
supplementary assumption leads to a further improvement 
of the fit of the mass-metallicity relation, but it is not an 
essential ingredient, because it alone cannot account for the 
mass-effective yield relation. 

Thus, the mai n cont ribution to the agreement with the 
ITremonti et al. I l|2004l ) data is the variable IGIMF, and not 
the details of the model. 

5.2 Star formation history 

We have chosen a rather simple recipe for the SFR and 
adopted a dependence of its time scale on galaxy mass 
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on the basis of matching the observed trend in the cur- 
rent gas fractions on galaxy mass. Chemical evolution mod- 
els with constant yield but gas infall remain rather close 
to the relation of the closed-box mo del, unless subjected 
to massive and rapid inflow of gas l|Koppen fc Edmunds! 
1 19991 : iKoppen fc Hensleil l2005l ) . Because oxygen is an ele- 
ment whose enrichment occurs with negligible time delays, 
its present abundance is predominantly determined by the 
gas fraction, but irrespective of the time dependence of the 
SFR. 

This applies also to our models with a variable IMF and 
continuous star formation: in massive galaxies with their 
high SFR the yield is rather insensitive to changes in the 
SFR (Fig. U). The instantaneous yield changes very little in 
the course of time, as evident e.g. in Fig. [7] only for low 
galaxy masses do the models with different a give different 
effective yields, because during its time evolution the low 
SFR results in a strong change of the yield. Although in 
the models the IGIMF depends both on time and galaxy 
mass, it is the dependence on the latter parameter which 
clearly dominates. Therefore one has a direct mapping of the 
SFR-yield relation into the mass-effective-yield relation, as 
the close similarity of the curves in Figs. 3] and [7] shows. 
Furthermore, models with a constant star formation time 
scale have the same gas fraction for all galaxy masses, hence 
the closed-box relation Z — — yln(/gas) provides a common 
scaling factor between yield and abundance (Fig. [6]). 

Therefore, as far as a continuous star formation is con- 
cerned, a comparison of its history predicted by our mod- 
els with observations does not constitute a critical test of 
their validity or of the variable IMF concept. The presence 
of burst-like forms of star formation and the consequences 
in conjunction with a variable IMF are addressed later in 
this section. But it worth noting that the time scale for 
star formation in the universe is found to be ab out 7 Gyr 
iGlazebrook et alllioO^ : iBrinchmann et al.lliooil ') and that 
models with such a value (i.e. a = 2) match quite well the ob- 
served mass-met allicity relation (Figs. [6] and [8} . The reason 
is the already mentioned simple mapping of the SFR-yield 
relation into the mass-met allicity relation: With a value of 
7 Gyr a closed-box model with linear SFR {x = 1) gives 
a present gas fraction of /g^s = 0.13. Since the gas frac- 
tion and thus the SFR decrease exponentially, a galaxy of 
Mgai = 1O^°M0 has an average SFR of Mgai(l-/gas)/13Gyr 
= 0.67 M0 /yr. From Fig.|4]one obtains a yield slightly above 
solar value. The closed-box relation then gives the oxygen 
abundance Z = — yln(/gas) as about twice the yield, as is 
seen in Fig. |6] 

Observational data are available to compare with our 
adopted relat ion between star for mation time scale and 
galaxy mass: iHeavens et al.l (120041 ) determine the stellar 
populations in nearly 100000 galaxies from the continua ob- 
served in the Sloan Digital Sky Survey and deduce the time 
dependence of the SFR in galaxies as a function of their 
mass. While galaxies with masses of less than lO^^M© have 
star formation increasing slightly with time, which might in- 
dicate of a very slow infall of gas, the SFR for e.g. IO^Mq 
initially is about 100 times higher and declines rapidly, with 
a time scale of about 1 ... 2 Gyr. This behaviour also implies 
that today's low mass galaxies remain richer in gas than the 
massive systems. Our recipe gives a similar value of 2.5 Gyr 
for IO^^Mq, but with 5 Gyr for lO^^M© the gas consumption 



is too rapid. For lower masses, our time scale increases and 
thus the star formation history is dominated by the infall 
which results in a nearly constant or increasi ng rate. This 
can be considered consistent with the data of iHeavens et al.l 
(|2004 ) which lump all masses bel ow 10^°Mq in a single bin. 

In a sample of the same size IBrinchmann et all (|2004D 
find that for masses below about lO^^Mo the SFR has a tight 
relation with the stellar mass in a galaxy (their Fig. 17). 
However, the high ratios of present and average past SFR 
confirm that star formation is dominated by bursts, and 
therefore its current value is not necessarily representative 
for the average rate of gas consumption. Therefore, for more 
detailed modelling of these galaxies the SFR recipe would 
need to include a proper description of the burst-type nature 
of star formation. 

We have refrained from doing this in the present work, 
but in Section [4.41 we illustrate the effects of starbursts on 
the results. In this extreme model with SFR/ < SFR > = 
30, which is about ten times larger than the largest value 
shown in Fig. 24 of IBrinchmann et al.l l|2004l ). the effective 
yield and the oxygen abundance for galaxies with < 10^" Mq 
are enhanced by about 0.4 dex. Since these galaxies have 
SFR < O.7M0 /yr, we see from Fig.Ujthat the yield increases 
by 0.3 ... 0.4 dex for an increase of a factor of 30 in the 
SFR, and 0.1 ... 0.2 dex for a factor of 3. Thus the enhance- 
ments would be less t han half of those s hown in Figs. 1131 
and 1141 The results of IBrinchmann et al. I (Hqm) 

would im- 
ply that while at masses near 10^° Mq there would a neg- 
ligible change from our results but that due to the effects 
of starbursts the effective yields and the oxygen abundances 
would progressively be enhanced for lower masses, amount- 
ing to about 0.2 dex at 1O*M0. The low-mass branch in both 
relations would be somewhat flatter than predicted from our 
continuous star formation models. However, this difference 
remains comparable with the uncertainties of the stellar nu- 
cleosynthesis prescriptions and of the observational data, to 
be discussed in the subsequent sections. 



5.3 Oxygen abundances 

The absolute values of abundances and yields from the 
models depend on the adopted recipes for the stellar nu- 
cleosynthesis, which may differ in the underlying stellar 
physics, such as the treatment of convection, the nuclear 
cross sections, and the explosion characteristics. For exam- 
ple, the stellar oxygen y i elds f rom Irhielemann et al.l (|l996f) 
and lWooslev fc Weaved ^99^) can differ for the same stel- 
lar mass by as much as a factor of two, eis shown by 
Thomas ct al. (1998i). Other effect s, such as stellar mass loss 
(|Maedei]|l992l ') and rotation (e.g. iMevnet fc Maederll2002t l 
also have a strong influence on the stellar yields. In view of 
these uncertainties as well as of those in the observational 
data, which are addressed below, the agreement of the abso- 
lute abundance value s between our models and the data of 
iTremonti et al.l (120041 ) should be regarded but as a welcome 
aspect. The essential feature of the variable IMF models is 
the agreement in the overall shape of both relations, already 
apparent in the closed box models. 
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5.4 Galaxy formation 



One basic assumption in the models is that the star forma- 
tion history of a galaxy is governed by an overall time scale. 
If galaxies form as isolated entities, one could expect such a 
behaviour. If - to take the hierarchical galaxy formation to a 
limiting case - all stars were formed in smaller substructures 
and in the same way, from which galaxies are built by mere 
additions, one should expect that all galaxies should have 
the same present gas fraction. The variation of gas fractions 
and present star formation among spiral and irregular galax- 
ies gives evidence that a substantial star formation must 
have taken place during or after the assembly of a galaxy. 
In the first place, a more massive galaxy would have under- 
gone a larger number of star formation episodes which could 
well result in an overall relation between star formation and 
present galactic mass. What fractions of st ars are formed in 
either of these two modes is not yet clear. IConselic3 (|2006l ) 
finds that most of today's massive galaxies probably have 
accreted the majority of their mass in only 4 or 5 major 
mergers in the early universe. As the merging objects already 
must have been quite massive (~ 2 x 10^" Mq in baryons) 
and young, the resulting IGIMF would be quite similar to 
the final IGIMF as the change in the IGIMF is most prom i- 
nent for galaxie s below 10^° Mm jWeidner fc KroupalliooH ) . 
lUnavane (liggeD and IVenn et all (|2004l ') have shown 

that the Milky Way stellar population has abundance pat- 
terns that cannot be derived from accreted satellites, unless 
they had significantly different stellar populations than to- 
day's dwarf spheroidal galaxies. This, however, would be 
unphysical, because 10^ Af© satellites would not be able 
to chemically enrich to the Milky Way values. 

If a major portion of the present stellar population was 
formed in smaller substructures, probably with a high SFR, 
or if each merging of a substructure results in an episode 
of enhanced star formation, one should expect higher val- 
ues for the average effective yield and the metallicity - as 
demonstrated in the models with bursty star formation his- 
tories. Thus the variation of the metallicities could be less 
than computed from our models. 

Our models - closed box and infall - thus constitute an 
upper limit of the effects due to a variation of the IGIMF 
on the SFR. 



5.5 Galactic Outflows 

The effective yields can further be reduced by the loss of 
metals due to galactic winds. While the details of the effects 
of supernova explosions on the chemical evolution are quit e 
complex (e.g. lMac Low fc Ferrar"alll999l : lRecchi et al.ll2006l ). 
one may distinguish two simplified cases: 

Some part of the metal-rich ejecta freshly produced by 
the massive stars are directly expelled from the galaxy. How 
large this fraction is depends on the details of the explosion, 
the multi-phase structure of the IS M, and the location in the 
galaxy (see e.g. iFrever et al.ll2003l ). which makes it difficult 
to give reliable estimates. It is obvious, however, that the 
effective yield could be reduced to arbitrarily small values, 
as it might be possible that all stellar ejecta leave the galaxy . 
Indee d, this issue is unclear, as lSihch fc Tenorio-TaglJ(|l998l . 
show that even massive galaxies would loose all their 
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Figure 15. Same as Fig. [8] but with additional gas outflows. 
The star formation time scale is 7.5 Gyr (a = 2), independent of 
galaxy mass. The top series of symbols (open ones) are models 
without outflows; for the ones below the ratio of outflow rate to 
SFR is taken to be constant at 1, 2, and 5. The symbols indicate 
the gas fractions, as in Fig. [6] 



metals through outflows if the iMac Low fc Ferraral (Il999l ) 
estimates are adopted. 

The other case is that the galaxy loses part of its ISM, 
following each event of star formation. Thus, the outflow 
rate would be expected to be proportional to the SFR. In 
Fig. [T5] we illustrate how the infall models of Fig. [S] are al- 
tered by outflows of this type: if the ratio of outflow and star 
formation rates is 1, 2, or 5, the metallicities are reduced by 
0.15, 0.3, and 0.75 dex, respectively. The effective yields are 
reduced by 0.3, 0.5, and 1.1 dex. Since in these models we 
use a constant star formation time scale, the gas fractions do 
not vary with galaxy mass. Thus, both the mass-metallicity 
and mass-yield relations are essentially modified by a mul- 
tiplicative factor. Note that the outfiows also will lower the 
gas fractions. 

The reduction factor of the effective yield depends 
strongly on the star formation time scale and the ratio of 
outfiow rate and SFR, but much less on the infall time scale, 
and the exponent x of the SFR, as presented in Fig. 1161 Since 
suitable models have star formation time scales of about 7 
Gyr (a = 2) , we may take from these schematic models that 
in order to cover the difference of about 1 dex in oxygen 
abundances and in the yields between massive and dwarf 
galaxies, the ratio of outfiow rate and SFR needs to vary 
only between and about 5. This number is likely to be 
near the lower end of this range, because stellar ejecta may 
also escape directly from dwarf galaxies. 

The combination of gas outflows and a variable IGIMF 
yields a mass-metallicity relation that is steeper than the 
one produced by either process alone. On the other hand, 
the presence of strongly fluctuating star formation in dwarf 
galaxies leads to higher yields - as shown in Section 14.41 
- which counteracts the effects of the variable IGIMF and 
thus gives a flatter mass-metallicity relation. Thus, chemi- 
cal evolution models with parametrised recipes would not be 
sufficient to discriminate between the three aspects. To fur- 
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Figure 16. The reduction factors of tlie effective yield below the 
true yield for models where the outflow rate is a constant multiple 
of the SFR. The upper set of curves refer to models with low SFR 
(a = 0.5), the lower one to those with a = 4. The full curves are 
without infall, the dashed curves are with infall time scales of 1, 
5, and 10 Gyr (shortest dashes). The dash-dotted curves are for 
an infall time scale of 5 Gyr, but with a quadratic SFR {x = 1). 

ther quantify this is beyond the scope of the present work, 
which focuses on the consequences of a variable IGIMF on 
the chemical properties of galaxies. 

5.6 Observational data 

To compare our models with the observed mass- 
metallicity relation w e use the recent observational data of 
iTremonti eral] l|2004l ) which constitutes a rather large and 
also homogeneous set of data. While our study is not the 
place to discuss in depth all uncertainties and limitations 
of this data set, several points need to be addressed: the 
abundances are estimated from the strong lines, thus with- 
out measuring the electron temperature in the ionised gas 
of each object which would require the observation of the 
faint diagnostic lines. While this can be used with calibration 
with photoionisation models to obtain oxygen abundances, 
direct comparison with determinations with measured elec- 
tron temperatures have shown that strong line methods may 
ove restimate the oxygen abundance by as much as 0.5 dex 
(cf. iKennicutt et al.ll200j ). Furthermore, the data pertain to 
the central regions of the galaxies which are likely to be of 
higher metallicity than the averaged value which the models 
are designed to compute. The computation of the effective 
yields requires the knowledge of the gas content of each ob- 
ject, which would need a determination of the mass in stars 
and gas in neutral, ionised, and molecular form. The deriva- 
tion of the gas fraction is further made more complicated 
in that such a measure should exclude any gas that does 
not take part in the star-gas cycle of evolution. Such details 
of the galaxies are no t available for the objects studied by 
ITremonti et al.l l|2004 ): they can determine the gas fraction 
but by the indirect means of determining the gas surface 
density from the Ha luminosit y via the relatio n between 
SFR and gas surface density by iKennicuttI l|l998l ). But this 



Figure 17. C omparison of th e stella r-mass-metallicity relation 
determined by ITremonti et al.l l|2004 ^ with a bundances for indi- 
vidua l gas-rich galaxies from the sample of iKoppen fc Henslerl 
l|2005h . Large circles indicate spiral galaxies, small filled circles 
are irregulars. The dash-dotted line is a fit to Garnett's (2002) 
data on oxygen abundances and H I mass. 

verifies that the obtained relation between galaxy mass and 
effective yields is in good agreement with the data from a 
number of individual galaxies available in t he literature. 

W e have compared the relations of ITremonti et al.l 
to data from individua l gas-rich galaxies, taken f rom 
the literature (for details see iKoppen fc Hensleil |2005| ) for 
different objects than used in their comparison. In most 
of these galaxies, abundances are determined by measur- 
ing the electron temperature of the H ll gas. As shown in 
Fig llTI these individual data show oxygen abundances sys- 
tematically lo wer by about 0.3 dex, which is in line with 
the findings of lKennicutt et al.l (|2003|) in M 101. T his is also 
in agreement with the results bv lYin et al] (12006^ who com- 
pared the oxygen abundances for 531 spectra of star-forming 
galaxies from the SDSS, based on the direct determina- 
tion of electron temperatures from the [O III] 5007/4363 
line ratio with the B ayesian estimates via strong lines by 
ITremonti et al.l l|2004l ). They find that in about half of the 
sample the strong-line technique overestimated the abun- 
dance by 0.3 dex (their Fig. 2), which is correlated to the 
N/0 abundance ratio. In the framework of our models, such 
an offset could be accounted for by a slightly higher value 
of /3 = 2.2, for instance. Such a change would still be well 
within the range of acceptable values for this parameter. 
Merely for compar ison, we also pl ot the relation of H I mass 
and metallicity bv lGarnettI (|2002l ). 

We also compare the effective oxygen y ields. Figure 1181 
shows that the mass-yield relation from ITremonti et al.l 
(|2004l ) is quite compatible with individual measurements, 
as was also shown by them by their choice of samples from 
the literature. A slight tendency to give hig her values is 
also apparent in fig. 8 of Trcmont i et al.l l|2004 ). There is an 
appreciable dispersion for a given galactic mass, especially 
for masses below about lO^M©. But despite this dispersion, 
there is a strong tendency for small galaxies to have also 
low effective yields. Such a feature may be the consequence 
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Figure 18. Co mparison of the stellar-mass— effective-yield rela- 
tion derived by iTremonti et alj ||2004|') with data for individua l 
gas-rich galaxies from the sample of iKoppen &: Hensler l|2005h. 



Large circles indicate the spiral galaxies in the sample of GarnettI 
||2002|). A filled sym bol indicates that the gas fraction is taken 
from lGarnet^ l|2002h . The dot-dashed curve is the mass— yield re- 
lation for the galactic wind model from ITremonti et al.l ||2004|) . 



Figure 19. The deviations of t he oxygen abun dances from the 
mean mass— metallicity relation o f iGarnettI l|2002il as a function of 
the deviations of the logarithmic gas fractions from the assumed 
relation shown in Fig.[9| Large c ircles indicate the spiral galaxies 
in the sample of iGarnetti | |2002|) . The dotted line indicates the 
slope of unity which would be expected for closed-box chemical 
evolution with the same constant effective yield. 



of a g alactic wind scenario as proposed by ITremonti et al.l 
l|2004 l as well as a signature for a variable IGIMF. One 
notes that most of the galaxies less massive than 10^ M© 
have higher effective yie lds than predict e d by the galactic 
wind model adopted by ITremonti et al.l (|2004 ) . The same 
holds true for our models, of course. Apart from larger un- 
certainties in the observational data expected from these 
dwarf galaxies, and the possibility that in these objects the 
true gas content could have been overestimated due to the 
presence of H I envelopes (cf. the example of Leo A dis- 
cussed bv lKoppen fc HensleiibOoH) . our model offers a pos- 
sible physical explanation for such a trend: bursts of intense 
star formation result in higher yields and higher metallicities 
than continuous star forming to the same present gas frac- 
tion. Nea rby dwarf gal axies have quite bursty star formation 
histories (|Grebelll200il l. which would tend to favour such an 
interpretation. Detailed studies of individual galaxies will be 
necessary to confirm this possibility. 

5.7 Scatter in the relations 

There is a substantial spread about the mass-metallicity re- 
lation, and an even larger one about the mass-gas-fraction 
relation. It is not the aim of this work to model this scat- 
ter, nor to provide physical explanations. But the chemical 
models permit to make the following statement: one could 
account for the spread in the gas fractions (Fig. llOp by in- 
voking, for instance, an intrinsic variation of a factor of 2 in 
the star formation time scale for galaxies of the same mass, 
as indicated by the three series of models. The spread of the 
oxygen abundances from these mode ls in Fig. | ll | is somewha t 
larger than the dispersion found bv ITremonti et all (|2004l ). 
Because of the relationship Z = yln(/gas) between metal- 
licity and fraction which exists in a closed-box model but is 
still closely valid in monotonic infall models, an increase in 



ln(/gas), caused by a decrease in the SFR, gives rise to an 
increase in metallicity. In our models with a variable IGIMF 
the lower SFR implies a smaller yield, thus the metallic- 
ity does not increase as much as if the yield were constant. 
While this can be seen in the results, it could constitute a 
possible test for the validity of our scenario: provided that 
gas fractions and abundances can be measured with suffi- 
cient accuracy so as to determine genuine differences in in- 
dividual objects, and if a correlation between deviations in 
abundances and gas fractions is found, our scenario would 
predict that the deviations in abundance should be less than 
in gas fractions. It is obvious that this quite subtle effect is 
not an easy test. 

If one assumed a relation between galactic H I mass and 
the gas fraction in our sample of individual galaxies, the de- 
viations from this relation are correlated with the deviations 
of the oxygon abundances from the mass-metallicity relation 
of Garnett (2002), as shown in Fig. 1191 This correlation is 
only marginally significant (correlation coefficient 0.6 for all 
objects and 0.5 for the spiral galaxies); but as the dotted 
line indicates, the trend is in good agreement to what one 
should expect if all galaxies had the same yield and evolved 
as nearly closed boxes. 



6 CONCLUSIONS 

The closed box model as well as the models with monoton- 
ically decreasing infall provide a remarkable match to the 
observed median mass-metallicity relation. This is in strong 
contrast to models with a constant IMF which all have the 
same effective yield for all galactic masses and all star for- 
mation time scales and therefore do not provide a reason- 
able match. In the case of the infall models also the effec- 
tive yields match the observed relation very reasonably. The 
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match with the observations is independent of the adopted 
infall time scale. For models with a starburst the result- 
ing IGIMF includes a higher number of massive stars. In 
comparison to continuous star formation the average metal 
yield will therefore be larger. This leads to to both higher 
oxygen abundances and higher effective yields in strong star- 
bursts. It must be noted that the rather schematic model 
of starbursts used here mainly shows the tendency of the 
results for extreme burst-type star formation histories. Re- 
gardless of the details it is important that these results pre- 
dict higher metallicities and higher oxygen yields for galaxies 
with strongly fluctuating star formation histories. Our mod- 
els - closed box and infall - thus constitute an upper limit of 
the effects due to a variation of the IGIMF on the SFR. By 
combining the effects of gas outflows and a variable IGIMF 
a steeper mass-met allicity relation would be achieved than 
the ones produced by either process alone. But in dwarf 
galaxies strong fluctuations in star formation would lead to 
higher yields thereby compensating this and flattening the 
stellar-mass-metallicity relation to the observed slope. 

It may be possible to test our model in more detail pro- 
vided that gas fractions and abundances can be measured 
with sufficient accuracy so as to determine genuine differ- 
ences in individual objects. If a correlation between devia- 
tions in abundances and gas fractions is found, our scenario 
would predict that the deviations in abundances should be 
less than in gas fractions. 

The idea of a SFR-dependent IGIMF provides a most 
attractive way to explain the observed relations between 
metallicity and galactic mass as well as the one between 
effective yield and galactic mass. The main effect at work 
is the lowering of the effective upper mass limit of stars for 
low star-formation rates, which reduces the number of su- 
pernovae of type II. 

Our results demonstrate a good agreement with the 
mass-metallicity data, even for (3 = 2, without needing to 
invoke outflows. This may suggest that the importance of 
outflows for shaping the mass-metallicity relation may need 
reconsideration. Thus, the oxygen yield becomes a sensitive 
function of the SFR but also of the cluster MF power-law 
index (5. Low-mass galaxies therefore appear chemically less- 
evolved than massive galaxies despite having the same ages. 

If one wanted to explain the mass-metallicity relation 
by an explicit variation of the slope of the IGIMF on galaxy 
mass, one would need to vary the exponent a above about 
1.5 Mq from 2.5 for galaxy masses of about 5 x 10* M© to 
2.3 for 2 X 10^^ Mq, which is the ran g e of t he data covered 
by the observations of lTremonti et all (|2004l ') . These changes 
of the slope IGIMF a correspond to t he 'minimal' scenario 
proposed in Weidner fc Kroupal (^^^. 



In a recent paper Lee et al] l|2006l ) investigate the mass- 
metallicity relation with a more recent sample of dwarf 
galaxies with Spitzer data. The sample is in principle in 
agreement with the data we use and form an extension to 
lower galactic masses. They flnd a rather flat mass-yield 
relation noting it "is difficult to explain, if galactic winds 
are ubiquitous in dwarf galaxies." Therewith they indicate 
that other solutions rather than galactic winds might still 
be possible. Here the flat stellar-mass-yield relation is ex- 
plainable with a SFR-dependent IGIMF by invoking bursty 

star-formation histories for dwarfs. 

Furthermore in another recent paper (|Dalcantonll2006l ') 



it has been shown that neither outflows nor accretion of 
low-metallicity gas can reproduce the low effective yields ob- 
served in low mass galaxies. It is conceivable that a variable 
IGIMF could help to explain this problem. 

Our models indicate that massive spirals and ellipti- 
cals should show less variations in the IGIMF and thus 
metallicities due to their high average SFR. Th is is ob- 
served in nearby galaxies as is seen in fig. 9 in iGarnettI 
(.2004) or in the samp le of 53000 galaxies plotted in fig. 6 in 
iTremonti l|2004) . The direct consequence of a depen- 

dence of the IGIMF on the SFR is a strong dependence of the 
oxygen abundance on the mass of a galaxy. It alone would 
already suffice to explain the observed mass-metallicity re- 
lation. Thus it would impose lesser demands on the loss of 
metals by outflows. 

To verify whether the SFR-dependence of the IGIMF 
could indeed be the origin or at least part of the explana- 
tion of the mass-metallicity relation of galaxies, we are in 
the process of comparing model predictions with other ob- 
servational properties, among them the abundances of other 
elements. The variation of the IGIMF affects strongly the 
massive stars but only weakly the intermediate mass stars. 
Thus we would expect that elements produced only in mas- 
sive stars will vary as strongly as oxygen, but that elements 
with a strong contribution from intermediate mass stars - 
such as carbon, nitrogen, and iron - will vary to a lesser 
extend. This will be the subject of a future study. 
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